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This is a paper about geometry and how one can derive several 
fundamental laws of physics from a simple postulate of geometrical 
nature. The method uses monogenic functions analysed in the alge- 
bra of 5 -dimensional spacetime, exploring the 4-dimensional waves 
that they generate. With this method one is able to arrive at equations 
of relativistic dynamics, quantum mechanics and electromagnetism. 
Fields as disparate as cosmology and particle physics will be influ- 
enced by this approach in a way that the paper only suggests. The pa- 
per provides an introduction to a formalism which shows prospects 
of one day leading to a theory of everything and suggests several 
areas of future development. 

1 Introduction 

The editor's invitation to write a chapter for this book about ether and the Uni- 
verse led me to think how my recent work had anything to do with ether, because 
the word was never used previously in my writings. It will become clear in the 
following sections that the concept of a privileged frame or absolute motion un- 
derlies all the argument. When one accepts the existence of a preferred frame, 
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the question of attaching that frame to some observable feature of the Universe is 
immediate. This question is addressed in Sec.[8]but we can anticipate that galaxy 
clusters are fixed and can be seen as the anchors for the preferred frame. This 
statement seems inconsistent with the observation that clusters of galaxies move 
relative to each other but it is resolved invoking an hyperspherical symmetry in 
the Universe that is revealed by the choice of appropriate coordinates. 

The relationship between geometry and physics is probably stronger in the 
General Theory of Relativity (GTR) than in any other physics field. It is the 
author's belief that a perfect theory will eventually be formulated, where geom- 
etry and physics become indistinguishable, so that the complete understanding 
of space properties, together with proper assignments between geometric and 
physical entities, will provide all necessary predictions, not only in relativistic 
dynamics but in physics as a whole. 

We don't have such perfect theory yet, however the author intends to show that 
GTR and Quantum Mechanics (QM) can be seen as originating from monogenic 
functions in the algebra of the 5-dimensional spacetime G41. These functions 
can generate a null displacement condition, thus reducing the dimensionality by 
one to the number of dimensions we are all used to. Besides generating GTR 
and QM, the same space generates also 4-dimensional Euclidean space where 
dynamics can be formulated and is quite often equivalent to the relativistic coun- 
terpart; Euclidean relativistic dynamics resembles Fermat's principle extended 
to 4 dimensions and is thus designated as 4-Dimensional Optics (4DO). 

Our goal is to show how the important equations of physics, such as relativity 
equations and equations of quantum mechanics, can be put under the umbrella 
of a common mathematical approachyj, |2D - We use geometric algebra as the 
framework but introduce monogenic functions with their null derivatives in order 
to advance the concept. Furthermore, we clarify some previous work in this 
direction and identify the steps to take in order to complete this ambitious project. 

Since A. Einstein formulated dynamics in 4-dimensional spacetime, this space 
is recognized by the vast majority of physicists as being the best for formulating 
the laws of physics. However, mathematical considerations lead to several alter- 
native 4D spaces. For example, the Euclidean 4-dimensional space of 4DO is 
equivalent to the 4D spacetime of GTR when the metric is static, and therefore 
the geodesies of one space can be mapped one-to-one with those of the other. 
Then one can choose to work in the space that is more suitable. We build upon 
previous work by ourselves and by other authors about null geodesies, regarding 
the condition that all material particles must follow null geodesies of 5D space: 

The implication of this for particles is clear: they should travel on 
null 5D geodesies. This idea has recently been taken up in the liter- 
ature, and has a considerable future. It means that what we perceive 
as massive particles in 4D are akin to photons in 5D.J3] 
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Accordingly, particles moving on null paths in 5D (dS 2 = 0) will 
appear as massive particles moving on timelike paths in 4D (ds 2 > 0) 

...it 



We actually improve on these null displacement ideas by introducing the more 
fundamental monogenic condition, deriving the former from the latter and estab- 
lishing a common first principle. 

The only postulates in this paper are of a geometrical nature and can be sum- 
marized in the definition of the space we are going to work with; this is the 4- 
dimensional null subspace of the 5-dimensional space with signature ( — h + + 
+). The choice of this geometric space does not imply any assumption for physi- 
cal space up to the point where geometric entities like coordinates and geodesies 
start being assigned to physical quantities like distances and trajectories. Some 
of those assignments will be made very soon in the exposition and will be kept 
consistently until the end in order to allow the reader some assessment of the 
proposed geometric model as a tool for the prediction of physical phenomena. 
Mapping between geometry and physics is facilitated if one chooses to work 
always with non-dimensional quantities; this is done with a suitable choice for 
standards of the fundamental units. From this point onwards all problems of 
dimensional homogeneity are avoided through the use of normalizing factors 
listed below for all units, defined with recourse to the fundamental constants: 
h — > Planck constant divided by 2n, G — > gravitational constant, c — > speed of 
light and e — > proton charge. 



This normalization defines a system of non-dimensional units (Planck units) 
with important consequences, namely: 1) All the fundamental constants, h, G, 
c, e, become unity; 2) a particle's Compton frequency, defined by v = mc 2 /h, 
becomes equal to the particle's mass; 3) the frequent term GMj (c 2 r) is simplified 
ioM/r. 

5-dimensional space can have amazing structure, providing countless parallels 
to the physical world; this paper is just a limited introductory look at such struc- 
ture and parallels. The exposition makes full use of an extraordinary and little 
known mathematical tool called geometric algebra (GA), a.k.a. Clifford algebra, 
which received an important thrust with the works of David Hestenes J5J. A 
good introduction to GA can be found in Gull et al. |@] and the following para- 
graphs use basically the notation and conventions therein. A complete course on 
physical applications of GA can be downloaded from the internet CD; the same 
authors published a more comprehensive version in book form 1 8]. An accessible 
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presentation of mechanics in GA formalism is provided by Hestenes 0]. This is 
the subject of first section, where some essential GA concepts and notation are 
introduced. 

Section two deals with monogenic function in flat 5D spacetime, deriving spe- 
cial relativity and the free particle Dirac equation from this simple concept. 4DO 
appears here as a perfect equivalent to special relativity, where trajectories can 
be understood as normals to 4-dimensional plane-like waves. The following sec- 
tion improves on this by allowing for curved space, introducing the notion of 
refractive index tensor. Section five examines the variational principle applied 
in both 4DO and GTR spaces to justify the equivalence of geodesies between 
the two spaces for static metrics. Refractive index is then related to its sources 
and the sources tensor is defined. The case of a central mass is examined and the 
links to Schwarzschild's metric are thoroughly discussed. Electromagnetism and 
electrodynamics are formulated as particular cases of refractive index in section 
seven and the sources tensor is here related to a current vector. The next sec- 
tion introduces the hypothesis of an hyperspherical symmetry in the Universe, 
which would call for the use of hyperspherical coordinates; the consequences 
for cosmology would include a complete dismissal of dark matter for a flat rate 
Hubble expansion. Before the conclusion, section nine shows how the mono- 
genic condition is effective in generating an SU (4) symmetry group and makes 
some advances towards a relation with the standard model of particle physics. 



2 Introduction to geometric algebra 

Geometric algebra is not usually taught in university courses and its presence in 
the literature is scarce; good reference works are 0.0,0]. We will concentrate 
on the algebra of 5-dimensional spacetime because this will be our main working 
space; this algebra incorporates as subalgebras those of the usual 3-dimensional 
Euclidean space, Euclidean 4-space and Minkowski spacetime. We begin with 
the simpler 5D flat space and progress to a 5D spacetime of general curvature. 

The geometric algebra G41 of the hyperbolic 5-dimensional space we consider 
is generated by the coordinate frame of orthonormal basis vectors o a such that 

(tf ) 2 = -l, 

(a,-) 2 = l, (2.1) 

<7 a -(7p=0, GC^jS. 

Note that the English characters i, j, k range from 1 to 4 while the Greek char- 
acters a, /3, 7 range from to 4. See the Appendix lAl for the complete notation 
convention used. 

Any two basis vectors can be multiplied, producing the new entity called a 
bivector. This bivector is the geometric product or, quite simply, the product; 
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this product is distributive. Similarly to the product of two basis vectors, the 
product of three different basis vectors produces a trivector and so forth up to the 
frvevector, because five is the dimension of space. 

We will simplify the notation for basis vector products using multiple indices, 
i.e. (7 a (Jp = <y a p- The algebra is 32-dimensional and is spanned by the basis 

• 1 scalar, 1, 

• 5 vectors, <7 a , 

• 10 bivectors (area), a aj g, 

• 10 trivectors (volume), o a py, 

• 5 tetravectors (4- volume), io a , 

• 1 pseudoscalar (5-volume), i = 001234- 
Several elements of this basis square to unity: 

(a ; ) 2 = (a 0l ) 2 = (o 0lJ ) 2 = (icr ) 2 = 1. (2.2) 

It is easy to verify the equations above; suppose we want to check that (oby) = 
1. Start by expanding the square and remove the compact notation (ob,j) 2 = 
OoOjOjOoOjOj, then swap the last Oj twice to bring it next to its homonymous; 
each swap changes the sign, so an even number of swaps preserves the sign: 
(oby) 2 = CTo<7i(<7/) 2 Oo<7;. From the third equation (12.11) we know that the squared 
vector is unity and we get successively (oo,j) 2 = OoOjOoOi = — (ob) 2 (a/) 2 = 
— (00) 2 ; using the first equation (12.11) we get finally (oo,j) 2 = 1 as desired. 

The remaining basis elements square to —1 as can be verified in a similar 
manner: 

(a ) 2 = {ajf = (o ljk f = (id) 2 = i 2 = -1. (2.3) 

Note that the pseudoscalar i commutes with all the other basis elements while 
being a square root of —1; this makes it a very special element which can play 
the role of the scalar imaginary in complex algebra. 

We can now address the geometric product of any two vectors a = a a o a and 
b = h^op making use of the distributive property 

ab= + + £ a a bVc a p; (2.4) 

V ' / cc^p 

and we notice it can be decomposed into a symmetric part, a scalar called the 
inner or interior product, and an anti-symmetric part, a bivector called the outer 
or exterior product. 

ab = a-b + aAb, ba = a-b — aAb. (2.5) 
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Reversing the definition one can write inner and outer products as 

a-b = —(ab + ba), a Ab = -(ab — ba). (2.6) 

The inner product is the same as the usual "dot product," the only difference 
being in the negative sign of the ao^o term; this is to be expected and is similar 
to what one finds in special relativity. The outer product represents an oriented 
area; in Euclidean 3-space it can be linked to the "cross product" by the relation 
cross(a,b) = —(7123a A b; here we introduced bold characters for 3-dimensional 
vectors and avoided defining a symbol for the cross product because we will not 
use it again. We also used the convention that interior and exterior products take 
precedence over geometric product in an expression. 

When a vector is operated with a multivector the inner product reduces the 
grade of each element by one unit and the outer product increases the grade 
by one. We will generalize the definition of inner and outer products below; 
under this generalized definition the inner product between a vector and a scalar 
produces a vector. Given a multivector a we refer to its grade-r part by writing 
< a > r ; the scalar or grade zero part is simply designated as < a >. By operating 
a vector with itself we obtain a scalar equal to the square of the vector's length 

a =aa = a-a-\-a/\a = a-a. (2.7) 

The definitions of inner and outer products can be extended to general multivec- 
tors 

a-b = J^(<a>a<b>p\ a _g,, (2.8) 

aAb = J E( <a> a <b> p)a+B- {2 - 9) 

Two other useful products are the scalar product, denoted as < ab > and commu- 
tator product, defined by 

axb = ab-ba. (2.10) 

In mixed product expressions we will always use the convention that inner and 
outer products take precedence over geometric products, thus reducing the num- 
ber of parenthesis. 

We will encounter exponentials with multivector exponents; two particular 
cases of exponentiation are specially important. If u is such that u 2 = — 1 and 
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is a scalar 



q2 q4 

= i-2, + 4 [ -...{=cose} + 

+u0-u— + ...{= usinO} (2.11) 



= cos0 + «sin0. 
Conversely if h is such that h 2 = 1 

= 1+M + | r+ »| r+ 9. + ... 
2 o 4 

+h0+h— + ...{= hsinhO} (2.12) 
= cosh0 +hsinhO. 

The exponential of bivectors is useful for defining rotations; a rotation of vector 
a by angle 6 on the On plane is performed by 

a ' = e°2iO/2 ae onO/2 = £ aR . (213) 

the tilde denotes reversion and reverses the order of all products. As a check we 
make a = 0\ 



e -a 12 fl/2 CTie ff 12 fl/2 = i C os - - on sin - I Oi * 



. e 

2" 2 

*(cos- + a 12 sin- (2.14) 

= cos0Ci + sin 0(72- 

Similarly, if we had made a = 02, the result would have been — sin 0o\ +cos 602- 
If we use B to represent a bivector whose plane is normal to Co and define its 
norm by \B\ = (BB) 1 / 2 , a general rotation in 4-space is represented by the rotor 

Rst ^ =m (my (2 , 5) 

The rotation angle is \B\ and the rotation plane is defined by B. A rotor is defined 
as a unitary even multivector (a multivector with even grade components only) 
which squares to unity; we are particularly interested in rotors with bivector 
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components. It is more general to define a rotation by a plane (bivector) then by 
an axis (vector) because the latter only works in 3D while the former is applicable 
in any dimension. When the plane of bivector B contains Oq, a similar operation 
does not produce a simple rotation but produces a boost, eventually combined 
with a rotation. Take for instance B = Oq\9/2 and define the transformation 
operator T = exp(5); a transformation of the basis vector Oq produces 

a' = To Q T = & -^ 6 l 1 a Q ^ 6 l 1 
= ^cosh ^ - Ooi sinh Ob * 

* ^cosh-^ + Oqi sinh-^ (2.16) 
= cosh 0(7q + sinh 9o\. 

In 5 -dimensional spacetime of general curvature, we introduce 5 coordinate 
frame vectors g a , the indices follow the conventions set forth in Appendix lAl 
We will also assume this spacetime to be a metric space whose metric tensor is 
given by 

gap =8a-gp; (2.17) 

the double index is used with g to denote the inner product of frame vectors and 
not their geometric product. The space signature is ( — (- + + +), which amounts 
to saying that goo < and ga > 0. A reciprocal frame is defined by the condition 



g a -gp = S a p. (2.18) 

Defining g a $ as the inverse of g a p, the matrix product of the two must be the 
identity matrix; using Einstein's summation convention this is 

g ar g Y p = 5 a p. (2.19) 
Using the definition (12.171) we have 

{g ar 8r)-gfi = S tt p; (2.20) 

comparing with Eq. (12.181) we determine g a with 

g a =g ay gy (2-21) 

If the coordinate frame vectors can be expressed as a linear combination of the 
orthonormed ones, we have 

ga=n^ a Op, (2.22) 

where n$ a is called the refractive index tensor or simply the refractive index, 
its 25 elements can vary from point to point as a function of the coordinates. |2j, 
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When the refractive index is the identity, we have g a = o a for the main or 
direct frame and g° = —(Jo, g l = <7, for the reciprocal frame, so that Eq. (12.181) is 
verified. In this work we will not consider spaces of general curvature but only 
those satisfying condition (12.221) . 

The first use we will make of the reciprocal frame is for the definition of two 
derivative operators. In flat space we define the vector derivative 

V = o a d a . (2.23) 

It will be convenient, sometimes, to use vector derivatives in subspaces of 5D 
space; these will be denoted by an upper index before the V and the particular 
index used determines the subspace to which the derivative applies; For instance 
"V = o m d m = o x d\ + a 2 ^2 + c 3 ^3. In 5-dimensional space it will be useful to 
split the vector derivative into its time and 4-dimensional parts 

V = -o d t + &di = -o d t + V. (2.24) 

The second derivative operator is the covariant derivative, sometimes called 
the Dirac operator, and it is defined in the reciprocal frame g a 

D = g a d a . (2.25) 

Taking into account the definition of the reciprocal frame (12.181) . we see that 
the covariant derivative is also a vector. In cases such as those we consider in 
this work, where there is a refractive index, it will be possible to define both 
derivatives in the same space. 

We define also second order differential operators, designated Laplacian and 
covariant Laplacian respectively, resulting from the inner product of one deriva- 
tive operator by itself. The square of a vector is always a scalar and the vector 
derivative is no exception, so the Laplacian is a scalar operator, which conse- 
quently acts separately in each component of a multivector. For 4+1 flat space 
it is 

V 2 = --^ + V 2 . (2.26) 
ot z 

One sees immediately that a 4-dimensional wave equation is obtained by zeroing 
the Laplacian of some function 

V 2 y/ = (- + V 2 ^) yr = 0. (2.27) 

This procedure will be used in the next section for the derivation of special rela- 
tivity and will be extended later to general curved spaces. 
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3 Monogenic functions and waves in flat space 



It turns out that there is a class of functions of great importance, called mono- 
genic functions, [8] characterized by having null vector derivative; a function y 
is monogenic if and only if 

Vy = 0. (3.1) 

A monogenic function has by necessity null Laplacian, as can be seen by dotting 
Eq. d3.il) with V on the left. We are then allowed to write 

Y d d li W = dmW- (3-2) 

i 

This can be recognized as a wave equation in the 4-dimensional space spanned 
by Oi which will accept plane wave type solutions of the general form 

¥ = ¥o ^ a + s ), (3.3) 

where yfQ is an amplitude whose characteristics we shall not discuss for now, 8 
is a phase angle and p a are constants such that 

I» 2 -(A)) 2 =0. (3.4) 

i 

By setting the argument of i// constant in Eq. (13.31) and differentiating we can 
get the differential equation 

p a dx a =0. (3.5) 

The first member can equivalently be written as the inner product of the two 
vectors p ■ dx = 0, where p = o a p a - In 5D hyperbolic space the inner product 
of two vectors can be null when the vectors are perpendicular but also when 
the two vectors are null; since we have established that p is a null vector, Eq. 
(13.51) can be satisfied either by dx normal to p or by (dx) 2 = 0. In the former 
case the condition describes a 3-volume called wavefront and in the latter case 
it describes the wave motion. Notice that the wavefronts are not surfaces but 
volumes, because we are working with 4-dimensional waves. 
The condition describing wave motion can be expanded as 

-(dx ) 2 + £(dx') 2 = 0. (3.6) 

This is a purely scalar equation and can be manipulated as such, which means 
we are allowed to rewrite it with any chosen terms in the second member; some 
of those manipulations are particularly significant. Suppose we decide to isolate 
(dx 4 ) 2 in the first member: (dx 4 ) 2 = (dx ) 2 — £(dx m ) 2 . We can then rename 
coordinate x 4 as r, to get the interval squared of special relativity for space-like 
displacements 

dr 2 = (dx ) 2 -£(dx m ) 2 . (3.7) 
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We have thus derived the space-like part of special relativity as a consequence 
of monogeneity in 5D hyperbolic space and simultaneously justified the physical 
interpretation for coordinates x° and x 4 as time and proper time, respectively. 

A different manipulation of Eq. (13.61) has great significance because it leads to 
the 4DO concept . 11 llll 1211 If we isolate (cbc ) 2 and replace x° by the letter t, we 
see that time becomes the interval in Euclidean 4D space 

df 2 = £(dx') 2 . (3.8) 

From this we conclude that the monogenic condition produces plane waves whose 
wavefronts are 3D volumes but can be represented by wavefront normals, just as 
it happens in standard optics with electromagnetic waves. 

Several readers may be worried with the fact that proper time is a line integral 
and not a coordinate in special relativity and so dr should not be allowed to ap- 
pear on the rhs of the equation. To this we will argue that the manipulations we 
have done, collapsing 5D spacetime into 4 dimensions through a null displace- 
ment condition and then promoting one of the coordinates into interval, is exactly 
equivalent to the process of defining a light cone in Minkowski spacetime and 
then applying Fermat's principle to define an Euclidean 3D metric on the light 
cone; we have just upgraded the procedure by including one extra dimension. 

The Dirac equation can also be derived from the monogenic condition but 
since it appears formulated in terms of matrices in all textbooks we will have 
to rewrite Eq. (13.11) also in terms of matrices, so that our GA manipulations can 
also be understood as matrix operations. This is easily achieved if we assign our 
frame vectors to Dirac matrices that square to the the identity matrix or minus 
the identity matrix as appropriate; the following list of assignments can be used 
but others would be equally effective 1 
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(3.9) 



There are 16 possible 4*4 Dirac matrices,! 13] of which we must choose 5 such that (<7o) 2 = 
— /, (t7j) 2 = / and cr a <7j3 = — GR<J a , for a^|5; the present choice will simplify our symmetry 
discussions further along. 
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There is no need to adopt different notations to refer to the frame vectors or to 
their matrix counterparts because the context will usually be sufficient to deter- 
mine what is meant. 

We can check that matrices o a form an orthonormal basis of 5D space by 
defining the inner product of square matrices as 

AB + BA 

AB= 1 . (3.10) 

It will then be possible to verify that the inner product of any two different o- 
matrices is null, (a ) 2 = — / and (cr') 2 = /; these are the conditions defining 
an orthonormal basis expressed in matrix form. A more formal approach to 
this subject would lead us to invoke the isomorphism between the complex al- 
gebra of 4 *4 matrices and Clifford algebra G4 1, the geometric algebra of 5D 



spacetime.[ 14]. 



It will now be convenient to expand the monogenic condition (13.11) as {o^d^ + 
c A d$) y/ = 0. If this is applied to the solution (13.31) and the derivative with respect 
to x 4 is evaluated we get 

(a% + aV)V = 0. (3.11) 

Let us now multiply both sides of the equation on the left by a 4 and note that 
matrix a 4 a° squares to the identity while the 3 matrices o A o m square to mi- 
nus identity; we rename these products as y-matrices in the form y" = o 4 o^. 
Rewriting the equation in this form we get 

(y%+ip 4 )V = 0. (3.12) 

The only thing this equation needs to be recognized as Dirac's is the replacement 
of p4 by the particle's mass m; simultaneously we assign the energy E to po and 
3D momentum p to <J m p m . 

We turn now our attention to the amplitude y/ in Eq. (13.31) because we know 
that the Dirac equation accepts solutions which are spinors and we want to find 
out their equivalents in our formulation. Applying the monogenic condition to 
Eq. ( 13.31) we see that the following equation must be verified 

Vo(° a Pa)=0. (3.13) 

If the as are interpreted as matrices, remembering that p is null, the only way the 
equation can be verified is by y/o being some constant multiplied by the matrix 
in parenthesis, which is a matrix representation of p. We can set the multiplying 
constant to unity and y/o becomes equal to p; the wavefunction y/ can then be 
interpreted as a Dirac spinor. The wave function in Eq. (13.31) can now be given a 
different form, taking in consideration the previous assignments 

y/ = A(o 4 m + p =f OQE y(±Et+p-x+mT+8) . (3 14) 
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where A is the amplitude and x = o m x m is the 3-dimensional position. 

In order to separate left and right spinor components we use a technique 
adapted from Ref. BSD. We choose an arbitrary 4x4 matrix which squares 
to identity, for instance 04, with which we form the two idempotent matrices 
(7+ 04) /2 and (/ — 04) /2. 1 These matrices are called idempotents because they 
reproduce themselves when squared. These idempotents absorb any 04 factor; 
as can be easily checked (1 + 04)04 = (I+04) and (1—04)04 = — (I— 04). 

Obviously we can decompose the wavefunction Xjf as 

I+04 I-04 
v = ¥ —— + v —— = v+ + v _. (3.15) 

This apparently trivial decomposition produces some surprising results due to 
the following relations 

e l6 (I+o 4 ) = (cos0+isin0)(7 + O4) 

= (7cose+ia 4 sin0)(7 + O4) (3.16) 
= e lc74e (/ + <7 4 ). 

and similarly 

e ie (/-a 4 )=e- ic74e (/-a4). (3.17) 

If we had chosen a different idempotent the result would have been similar; we 
will see how the various idempotents are arranged in a symmetry group and it 
has been argued that they may be related to elementary particles. 111511 



4 Relativistic dynamics 

When working in curved spaces the monogenic condition is naturally modified, 
replacing the vector derivative V with the co variant derivative D. A generalized 
monogenic function is then a function that verifies the equation 

D\jf = 0. (4.1) 

Similarly to what happens in flat space, the covariant Laplacian is a scalar and a 
monogenic function must verify the second order differential equation 

DV = 0. (4.2) 

It is possible to write a general expression for the covariant Laplacian in terms of 
the metric tensor components (see Section 2.1 1]) but we will consider only 
situations where that complete general expression is not needed. 

Matrix a is the same as matrix y 5 = iy°y l y 2 y i . 
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When Eq. ( 14. II) is multiplied on the left by D, we are applying second deriva- 
tives to the function, but we are simultaneously applying first order derivatives to 
the reciprocal frame vectors present in the definition of D itself. We can simplify 
the calculations if the variations of the frame vectors are taken to be much slower 
than those of function i// so that frame vector derivatives can be neglected. With 
this approximation, the covariant Laplacian becomes D 2 = g a ^d a p and Eq. (14.21) 
can be written 

8 aP dapY = 0. (4.3) 

This equation can have a solution of the type given by Eq. (13.31) if again the 
derivatives of p a are neglected. This approximation is usually of the same order 
as the former one and should not be seen as a second restriction. Inserting Eq. 
(13.31) one sees that it is a solution if 

g al5 PaPp=0. (4.4) 

This equation is the curved space equivalent to Eq. (13.41) and it means that the 
square of vector p = g a p a is zero, that is, p is a vector of zero length; for this 
reason it is called a null vector or nilpotent. Vector p is the momentum vector 
and should not be confused with 4-dimensional conjugate momentum vectors 
defined below. 

We arrive again at Eq. (13.51 ) and the condition describing 4D wave motion can 
be expanded as 

g a pdx a dxl 3 =0. (4.5) 

This condition effectively reduces the spatial dimension to four but the resulting 
space is non-metric because all displacements have zero length. We will remove 
this difficulty by considering two special cases. First let us assume that vector go 
is normal to the other frame vectors so that all go; factors are zeroed; condition 
(14.51) becomes 

g00 (dx ) 2 + gij dx l dxJ =0. (4.6) 

All the terms in this equation are scalars and we are allowed to rewrite it with 
(dx ) 2 in the lhs 

(cU ) 2 = -^dx'dx j . (4.7) 
goo 

We could have arrived at the same result by defining a 4-dimensional displace- 
ment vector 

dx°v=^Lg l dx i ; (4.8) 
v^oo 

and then squaring it to evaluate its length; v is a unit vector called velocity be- 
cause its definition is similar to the usual definition of 3-dimensional velocity; 
its components are 

dx' 
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Being unitary, the velocity can be obtained by a rotation of the 04 frame vector 



v = Ro 4 R. (4.10) 

The rotation angle is a measure of the 3-dimensional velocity component. A 
null angle corresponds to v directed along 04 and null 3D component, while a 
n/2 angle corresponds to the maximum possible 3D component. The idea that 
physical velocity can be seen as the 3D component of a unitary 4D vector has 
been explored in several papers but see [ 11711 . 

Equation ( 14.81) projects the original 5-dimensional space into an Euclidean 
signature 4 dimensional space, where an elementary displacement is given by the 
variation of coordinate x°. In the particular case where go = Oq the displacement 
vector simplifies to dx°v = gidx 1 and we can see clearly that the signature is 
Euclidean because the four g ; have positive norm. Although it has not been 
mentioned, we have assumed that none of the frame vectors is a function of 
coordinate x°. 

Returning to Eq. (14.61) we can now impose the condition that g4 is normal to 
the other frame vectors in order to isolate (dx 4 ) 2 instead of (dx ) 2 , as we did 
before; 

(dx*) 2 = -^dx»dx V . (4n) 
£44 

We have now projected onto 4-dimensional space with signature (H ), 

known as Minkowski signature. In order to check this consider again the special 
case with go = Ob an d the equation becomes 

(d /)2 = _L (d* ) 2 - ^dxW; (4.12) 

£44 g44 

the diagonal elements gu are necessarily positive, which allows a verification 
of Minkowski signature. Contrary to what happened in the previous case, we 
cannot now obtain (dx 4 ) 2 by squaring a vector but we can do it by consideration 
of the bivector 

dx 4 v= l_ g/1 gV. (4.13) 

All the products g^g 4 are bivectors because we imposed g4 to be normal to the 
other frame vectors. When (dx 4 ) 2 is evaluated by an inner product we notice that 
gog 4 has positive square while the three g m g 4 have negative square, ensuring that 
a Minkowski signature is obtained. Naturally we have to impose the condition 
that none of the frame vectors depends on x 4 . Bivector v is such that v 2 = vv = 1 
and it can be obtained by a Lorentz transformation of bivector Oo4- 

v = fa 04 r, (4.14) 
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where T is of the form T = exp(5) and B is a bivector whose plane is normal to 
04. Note that T is a pure rotation when the bivector plane is normal to both Oq 
and 04. 

In special relativity it is usual to work in a space spanned by an orthonormed 
frame of vectors 7^ such that (y ) 2 = 1 and (fm) 2 = — 1, producing the desired 
Minkowski signature [8]. The geometric algebra of this space is isomorphic to 
the even sub-algebra of G4J and so the area element dx 4 v (14.131) can be reformu- 
lated as a vector called relativistic 4- velocity. The four 7 bivectors are defined in 
a similar way to the 7 matrices used in Eq. (13.121) . which is to be expected from 
the isomorphism between geometric and matrix algebras already mentioned. 

Equations (14.71) and (14.1 II ) define two alternative 4-dimensional spaces, those 
of 4 -dimensional optics (4DO), with metric tensor —gij/goo and general theory 
of relativity (GTR) with metric tensor — guv/gAA, respectively; in the former x° 
is an affine parameter while in the latter it is x 4 that takes such role. In fact 
Eq. (14.1 II) only covers the spacelike part of GTR space, because (cbc 4 ) 2 is nec- 
essarily non-negative. Naturally there is the limitation that the frame vectors 
are independent of both x° and x 4 , equivalent to imposing a static metric, and 
also that goi = g/i4 = 0- Provided the metric is static, the geodesies of 4DO can 
be mapped one-to-one with spacelike geodesies of GTR and we can choose to 
work on the space that best suits us for free fall dynamics. For a physical in- 
terpretation of geometric relations it will frequently be convenient to assign new 
designations to the 5D coordinates that acquire the role of affine parameter in 
the null subspace. We recall the assignments x° = t and x 4 = v, total derivatives 
with respect to these coordinates will receive a special notation: df/dt = f and 
df/dz = f. 

Unless otherwise specified, we will assume that the frame vector associated 
with coordinate x° is unitary and normal to all the others, that is go = Co and 
go/ = 0. Recalling from Eq. (14.71) . these conditions allow the definition of 4DO 
space with metric tensor gjj. Although we could try a more general approach, 
we would loose the possibility of interpreting time as a line element and this, 
as we shall see, provides very interesting and novel interpretations of physics' 
equations. In many cases it is also true that #4 is normal to the other frame vectors 
and we have seen that in those cases we can make metric conversions between 
GTR and 4DO; as we shall see, electromagnetism requires a non-normal #4 and 
so we leave this possibility open. 

For the moment we will concentrate on isotropic space, characterized by or- 
thogonal refractive index vectors g, whose norm can change with coordinates 
but is the same for all vectors. Normally we relax this condition by accepting 
that the three g m must have equal norm but g4 can be different. The reason for 
this relaxed isotropy is found in the parallel we make with physics by assigning 
dimensions 1 to 3 to physical space. Isotropy in a physical sense need only be 
concerned with these dimensions and ignores what happens with dimension 4. 
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We will therefore characterize an isotropic space by the refractive index frame 
go = Ob, 8m = n r o m , g4 = W4O4. Indeed we could also accept a non-orthogonal 
g4 within the relaxed isotropy concept but we will not do so for the moment. 
Equation (14.71) can now be written in terms of the isotropic refractive indices 

as 

df 2 = (/v) 2 £(dx m ) 2 + (n 4 dT) 2 . (4.15) 

m 

Spherically symmetric static metrics play a special role; this means that the re- 
fractive index can be expressed as functions of r if we adopt spherical coordi- 
nates. The previous equation then becomes 

- - (n r ) 2 [dr 2 + r 2 (d0 2 + sin 2 0d<p 2 )] + (n 4 dx) 2 . (4.16) 



dt z 



Since we have g4 normal to the other vectors we can apply metric conversion 
and write the equivalent quadratic form for GTR 

2 2 

dT 2 =(^£) -(^j [dr 2 + r 2 (de 2 + sin 2 0d<p 2 )]. (4.17) 

In the case of a central mass, we can examine how the Schwarzschild metric 
in GTR can be transposed to 4DO. The usual form of the metric is 

-l 



dx 2 



1 



2M\ 
X ) 



I di- i 1 

2 a Am 2 



2M\ 

YJ 



dX 2 - 



-X 2 (d0 2 + sin 2 0d<p 2 ) ; 



(4.18) 



where M is the spherical mass and % is the radial coordinate, not the distance to 
the centre of the mass. This form is non-isotropic but a change of coordinates 
can be made that returns the expression to isotropic form (see DTnverno IliL 
section 14.7]): 



(x-M+^/x 2 -2M X S j /2; 



(4.19) 



and the new form of the metric is 

2 



(\ 



dx 1 



M ^ 

27 



M 



dt 2 



\ 1 + 2rJ 



M 



1 + — ) [dr 2 - r 2 (dd 2 + sin 2 0d<p 2 ) ] . (4.20) 



From this equation we immediately define two coefficients, which are called 
refractive index coefficients, 



«4 



M ' 



1 + 



M 
~2r 



1- 



2r 



M 
~2r 



(4.21) 
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These refractive indices provide a 4DO Euclidean space equivalent to Schwarzschild 
metric, allowing 4DO to be used as an alternative to GTR. Recalling that we de- 
rived trajectories from solutions (13.31) of a 4-dimensional wave equation (14.31) . 
it becomes clear that orbits can also be seen as 4-dimensional guided waves by 
what could be described as a 4-dimensional optical fibre. Modes are to be ex- 
pected in these waveguides and we shall say something about them later on. 

5 Fermat's principle in 4 dimensions 

Fermat's principle applies to optics and states that the path followed by a light 
ray is the one that makes the travel time an extremum; usually it is the path that 
minimizes the time but in some cases a ray can follow a path of maximum or 
stationary time. These solutions are usually unstable, so one takes the view that 
light must follow the quickest path. In Eq. (14.71) we have defined a time interval 
associated with a 4-dimensional elementary displacement, which allows us to 
determine, by integration, a travel time associated with displacements of any 
size along a given 4-dimensional path. We can then extend Fermat's principle 
to 4D and impose an extremum requirement in order to select a privileged path 
between any two 4D points. Taking the square root to Eq. (14.71) 

d/ = l-ilLfaifaj. (5.1) 
V £oo 

Integrating between two points Pi and P% 

t= J- — dx'dxJ = [ P \l-— xWdi. (5.2) 

JP\ V ^00 JPx V £00 

In order to evaluate the previous integral one must know the particular path link- 
ing the points by defining functions x l (t), allowing the replacement dx l = x'dt. 
At this stage it is useful to define a Lagrangian 

L=-^-x i x i . (5.3) 
2goo 

The time integral can then be written 

rPi , 

t= / V2Ldt. (5.4) 
JPx 

Time has to remain stationary against any small change of path; therefore 
we envisage a slightly distorted path defined by functions x l (t) +ex'(t), where 
£ is arbitrarily small and % l {t) are functions that specify distortion. Since the 
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distortion must not affect the end points, the distortion functions must vanish at 
those points. The time integral will now be a function of £ and we require that 



dt(e) 
de 



0. 



(5.5) 



£=0 



Now, the Lagrangian (15.31) is a function of x', through g a p and also an explicit 
function of x l . Allowing for a path change, through £ makes t in Eq. (15.41) a 
function of £ 



f(e) 



Pi 



P\ 



2L(x* + ex 1 + x l + £%>) dt. 



This can now be derived with respect to £ 



<fc(e) 
d£ 



e=0 



Pi 1 



Pi 



dL 



dL 



(5.6) 



(5.7) 



e=0 



Note that the first term on the rhs can be written 



Pl 1 dL 'iA 



VlLdx 
This can be integrated by parts 



p 2d(V2L) . iA 
\ dx' 



(5.8) 



*d(V2L) , 



Pi 



dx 1 



d(V2L) , 

* - A- 



ft 



p 2 d_ / ^(V2L ) 



Z'dr. (5.9) 



The first term on the second member is zero because X' vanishes for the end 
points; replacing in Eq. ( 15.71) 



dt(e) 


1 r p 2 


\-( 


1 dL\ 


d£ 


£=0 V2 Jpi 


dt V 


VZdx 1 ) 



+ 



1 dL 

VLdx 1 



X l dt. 



(5.10) 



The rhs must be zero for arbitrary distortion functions % l , so we conclude that 
the following set of four simultaneous equations must be verified 



1 dL 



1 dL 



dtyy/Ldx'J y/Ldx?' 



(5.11) 



these are called the Euler-Lagrange equations. 

Consideration of Eqs. (14.81) and (14.1 II) allows us to conclude that the La- 
grangian defined by (15.31) can also be written as L = v 2 /2 and must always equal 
1/2. From the Lagrangian one defines immediately the conjugate momenta 



dx 1 



X J . 



goo 



(5.12) 
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Notice the use of the lower index (v ; ) to represent momenta while velocity com- 
ponents have an upper index (v ! ). The conjugate momenta are the components 
of the conjugate momentum vector 

s'Vi 

(5.13) 



V-£oo 
and from Eq. (12.181) 



V-goov = g'vi = g'gijX J = g jx J . (5.14) 

The conjugate momentum and velocity are the same but their components are 
referred to the reciprocal and refractive index frames, respectively. 1 Notice also 
that by virtue of Eq. (13.41) it is also 

Vi = ^. (5.15) 
Po 

The Euler-Lagrange equations (15.1 II) can now be given a simpler form 

Vj = djL. (5.16) 

This set of four equations defines trajectories of minimum time in 4DO space as 
long as the frame vectors g a are known everywhere, independently of the fact 
that they may or may not be referred to the orthonormed frame via a refractive 
index. By definition these trajectories are the geodesies of 4DO space, spanned 
by frame vectors g,-/ yj— goo, with metric tensor —gij/goo- 

Following an exactly similar procedure we can find trajectories which extrem- 
ize proper time, defined by taking the positive square root of Eq. (14.1 II ). The 
Lagrangian is now defined by 

j? = -li»l^x v . (5.17) 

2 g44 



Consequently the conjugate momenta are 

v M = || = ^r. (5.18) 
dxV g u 



From Eq. (13.41) we have = Pfi/p4', the associated Euler-Lagrange equations 
are 

v M = ^. (5.19) 

"These are, by definition, spacelike geodesies of GTR with metric tensor — g^ v /g44 
and we have thus defined a method for one-to-one geodesic mapping between 
4DO and spacelike GTR. Recalling the conditions for this mapping to be valid, 
all the frame vectors must be independent of both t and x and go and g4 must 
be normal to the other 3 frame vectors. In tensor terms, all the g a p must be 
independent from t and x and go/ = gfiA — 0." 



In most cases goo — — 1 , tne velocity can be conveniently written v = gjX? and conjugate momenta 

Vi = gijX j . 
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6 The sources of refractive index 



The set of 4 equations (15.161) defines the geodesies of 4DO space; particularly in 
cases where there is a refractive index, it defines trajectories of minimum time 
but does not tell us anything about what produces the refractive index in the 
first place. Similarly the set of equations ( 15.191 ) defines the geodesies of GTR 
space without telling us what shapes space. In order to analyse this question we 
must return to the general case of a refractive frame g a without other impositions 
besides the existence of a refractive index. 
Considering the momentum vector 

P = Pag a = Panp V, (6.1) 

with riaJnP y = 8a , we will now take its time derivative. Using Eq. (IB .41) 

p = x-{Dp) = x-G (6.2) 

By a suitable choice of coordinates we can always have g° = o°. We can then 
invoke the fact that for an elementary particle in flat space the momentum vector 
components can be associated with the concepts of energy, 3D momentum and 
rest mass as p = Eo° + p + mo 4 (see Sec.|3]) If this consequence is extended to 
curved space and to mass distributions, we write p = Eo° + p + mg 4 , where now 
E is energy density, p = p m g m is 3D momentum density and m is mass density. 
The previous equation then becomes 

Eo°+p + mg 4 = x-G. (6.3) 

When the Laplacian is applied to the momentum vector the result is still nec- 
essarily a vector 

D 2 p = S. (6.4) 
Vector S is called the sources vector and can be expanded into 25 terms as 

S = (D 2 nP a )o p p a = sP a o pP a ; (6.5) 

where p a = g a P pp . Tensor S a p contains the coefficients of the sources vector 
and we call it the sources tensor. The sources tensor influences the shape of 
geodesies as we shall see in one particularly important situation. One important 
consequence that we don't pursue here is that by zeroing the sources vector one 
obtains the wave equation D 2 p = 0, which accepts gravitational wave solutions. 

If c° is normal to the other frame vectors we can write p = E(o° + v) in the 
reciprocal frame, with v a unit vector or p = E(—Oq + v) in the direct frame. 
Equation ( 16.21) can then be given the form 

E(o° + v) +Ev = Oq + v-G. (6.6) 
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Since G can have scalar and bivector components, the scalar part must be re- 
sponsible for the energy change, while the bivector part rotates the velocity v. 
The bivector part of G is generated by D A p, which allows a simplification of 
the previous equation to 

v = v-(DAv), (6.7) 

if the frame vectors are independent of t. This equation is exactly equivalent to 
the set of Euler-Lagrange equations ( 15.161 ) but it was derived in a way which tells 
us when to expect geodesic movement or free fall. 

We will now investigate spherically symmetric solutions in isotropic condi- 
tions defined by Eq. (14.161) : this means that the refractive index can be expressed 
as functions of r. The vector derivative in spherical coordinates is of course 

D = — (o r d r + - o e d e h )— o<pdm) - o t d t + — o T d T . (6.8) 

n r \ r rsmtJ J n$ 

The Laplacian is the inner product of D with itself but the frame vectors' deriva- 
tives must be considered; all the derivatives with respect to t, r and x are zero 
and the non-zero ones are 

deO,=Oe, d<pO r = sinflcfy, 

d e O0 = -<j r , dyCe = cos0a<p, (6.9) 
de = 0, dyOy = — sin 9o r — cos 9 Oq . 

After evaluation the curved Laplacian becomes 

D = --2 [d rr + -d r -^d r +^d ee + 
\n r ) V r n r ' 



COt0 csc 2 \ 1 
H — ~2~ oq H ^— dyy \-dn + j—y: vtt- (6.10) 

The search for solutions of Eq. (16.41) must necessarily start with vanishing 
second member, a zero sources situation, which one would implicitly assign to 
vacuum; this is a wrong assumption as we will show. Zeroing the second member 
implies that the Laplacian of both n r and must be zero; considering that they 
are functions of r we get the following equation for n r 

„ In' (n') 2 
+ — -^- = 0, (6.11) 
r n r 

with general solution n r = bexp(a/r). It is legitimate to make b = 1 because the 
refractive index must be unity at infinity. Using this solution in Eq. ( 16.101 ) the 
Laplacian becomes 



D 2 = e- a / r (d rr + -d r + ^d r + \d ee + 
\ r r L r L 



cotfl . csc 2 . \ _ 1 _ 
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which produces the solution H4 = n r . So space must be truly isotropic and not 
relaxed isotropic as we had allowed. The solution we have found for the refrac- 
tive index components in isotropic space can correctly model Newton dynamics, 
which led the author to adhere to it for some time W17I1 . However if inserted into 
Eq. (14.1 II) this solution produces a GTR metric which is verifiably in disagree- 
ment with observations; consequently it has purely geometric significance. 

The inadequacy of the isotropic solution found above for relativistic predic- 
tions deserves some thought, so that we can search for solutions guided by the 
results that are expected to have physical significance. In the physical world we 
are never in a situation of zero sources because the shape of space or the ex- 
istence of a refractive index must always be tested with a test particle. A test 
particle is an abstraction corresponding to a point mass considered so small as 
to have no influence on the shape of space; in reality a point particle is a black 
hole in GTR, although this fact is always overlooked; one wonders how a black 
hole is postulated not to influence space geometry. A test particle must be seen 
as source of refractive index itself and its influence on the shape of space should 
not be neglected in any circumstances. If this is the case the solutions for van- 
ishing sources vector may have only geometric meaning, with no connection to 
physical reality. 

The question is then what should go into the second member of Eq. (16.41) 
in order to find physically meaningful solutions. If we are testing gravity we 
must assume some mass density to suffer gravitational influence; this is what is 
usually designated as non-interacting dust, meaning that some continuous distri- 
bution of non-interacting particles follows the geodesies of space. Mass density 
is expected to be associated with 5 4 4; on the other hand we are assuming that 
this mass density is very small and so we use flat space Laplacian to evaluate it. 
We consequently make an ad hoc proposal for the sources vector in the second 
member of Eq. (16.41 ) 

S=-V 2 n 4 04- (6.13) 

Equation (16.41) becomes 

D 2 i = -V 2 n 4 a 4 ; (6.14) 
as a result the equation for n r remains unchanged but the equation for n 4 becomes 

" . 2n 4 n' r n' 4 „ 2n' 4 

n 4 -\ - = — n 4 -\ -. (6.15) 

r n r r 

When n r is given the exponential form found above, the solution is n$ = ^fn~ r . 
This can now be entered into Eq. (14.1 II) and the coefficients can be expanded 
in series and compared to Schwarzschild's for the determination of parameter a. 
The final solution, for a stationary mass M is 

n,. = e 2M /'', n 4 = e M/ '". (6.16) 
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The equivalent GTR space is characterized by the quadratic form 

dr 2 = e - 2M / r dt 2 -e 2M /'X((k'*) 2 . (6.17) 

m 

Expanding in series of M/r the coefficients of this metric one would find that 
the lower order terms are exactly the same as for Schwarzschild's and so the 
predictions of the metrics are indistinguishable for small values of the expansion 
variable. Montanus arrives at the same solutions with a different reasoning; 
Yilmaz was probably the first author to propose this metric lEoLEliE^ . 

Equation (16.141 ) can be interpreted in physical terms as containing the essence 
of gravitation. When solved for spherically symmetric solutions, as we have 
done, the first member provides the definition of a stationary gravitational mass 
as the factor M appearing in the exponent and the second member defines iner- 
tial mass as V 2 n^. Gravitational mass is defined with recourse to some particle 
which undergoes gravitational influence and is animated with velocity v and in- 
ertial mass cannot be defined without some field acting upon it. Complete 
investigation of the sources tensor elements and their relation to physical quan- 
tities is not yet done; it is believed that 16 terms of this tensor have strong links 
with homologous elements of stress tensor in GTR, while the others are related 
to electromagnetic field. 



7 Electromagnetism in 5D spacetime 

Maxwell's equations can easily be written in the form of Eq. (16.41) if we don't 
impose the condition that #4 should remain normal the other frame vectors; as 
we have seen in section this has the consequence that there will be no GTR 
equivalent to the equations formulated in 4DO. 

We will consider the non-orthonormed reciprocal frame defined by 

= < g 4 = ^-A^ + a 4 ; (7.1) 
m 

where q and m are charge and mass densities, respectively, and A =A fl o^ L is the 
electromagnetic vector potential, assumed to be a function of coordinates t and 
x m but independent of x. The associated direct frame has vectors 

q 

gu = On A„04, g4 = 04; (7.2) 

m 

and one can easily verify that Eq. (12.181 ) is obeyed. The momentum vector in 
the reciprocal frame is p = Eo° + p m o m + qA^o* 1 + mo 4 and G in the second 
member of Eq. ( 16.21) is G = qDA. We will assume D A to be zero, as one 
usually does in electromagnetism; also D can be replaced by because the 
vector potential does not depend on x. It is convenient to define the Faraday 
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bivector F = ^VA, similarly to what is done in Ref. the dynamics equation 
then becomes 

p + qA = qx-F; (7.3) 

and rearranging 

p = qx-F-qA. (7.4) 

The first term in the second member is the Lorentz force and the second term is 
due to the radiation of an accelerated charge. 

Recalling the wave displacement vector Eq. (IB.ll) we have now 

dx = o a dx a --A u o 4 dx^. (7.5) 
m 

This corresponds to a refractive index tensor whose non-zero terms are 

n a a = h n\ = -^-A ll . (7.6) 
m 

According to Eq. (16.51) the sources tensor has all terms null except for the 
following 

S% = -^D 2 A M ; (7.7) 
m 

where D is the covariant derivative given by 

D = g a d a = a% + (a 4 + ^-A^)d 4 . (7.8) 

m 

We can then define the current vector 7 verifying 

M V 2 a = ^VF = 7, (7.9) 

where 

7=--SV M - (7.10) 

q 

Please refer to |8, Chap. 7] or to |0, Part 2] to see how these equations generate 
classical electromagnetism. 

In free space we make 7 = and Eq. ( 17.91) accepts plane wave solutions for 
F which are of course electromagnetic waves. Notice that these solutions prop- 
agate in directions normal to proper time, which is perfectly consistent with the 
classical relativistic formulation. 

The Dirac equation for a free particle has been derived from the 5-dimensional 
monogenic condition in Sec.|3]but we are now in position to include the effects 
of an EM field. Because we are working in geometric algebra, our quantum 
mechanics equations will inherit that character but the isomorphism between the 
geometric algebra of 5D spacetime, G4 1, and complex algebra of 4*4 matrices, 
M(4,C), ensures that they can be translated into the more usual Dirac matrix 
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formalism. Electrodynamics can now be implemented in the the same way used 
in Sec. |7] to implement classical electromagnetism. The monogenic condition 
must now be established with the covariant derivative given by Eq. (17.81) 

<7%v/+ (a 4 + -A fl a^d 4 \if = 0. (7.11) 

Multiplying on the left by a 4 and taking d^y = \m\\f 

[f(d^+iqA fl )+im]\i/=0. (7.12) 

This equation can be compared to what is found in any quantum mechanics text- 
book.. 

It is now adequate to say a few words about quantization, which is inherent 
to 5D monogenic functions. We have already seen that these functions are 4- 
dimensional waves, that is, they have 3-dimensional wavefronts normal to the 
direction of propagation. Whenever the refractive index distribution traps one of 
these waves a 4-dimensional waveguide is produced, which has its own allowed 
propagating modes. In the particular case of a central potential, be it an atom's 
or a galaxy's nucleus, we expect spherical harmonic modes, which produce the 
well known electron orbitals in the atom and have unknown manifestations in a 
galaxy. 

8 Hyperspherical coordinates 

Deriving physical equations and predictions from purely geometrical equations 
is an exercise whose success depends on the correct assignment of coordinates 
to physical entities; the same space will produce different predictions if differ- 
ent options are taken for coordinate assignment. In the previous sections we 
assumed that empty space could be modelled by an assignment of time, three 
spatial directions and proper time to five orthogonal directions in 5D spacetime. 
We are now going to experiment with a different assignment of flat space coor- 
dinates, which will explore the possibility that physics and the Universe have an 
inbuilt hyperspherical symmetry. The exercise consists on assigning coordinate 
x 4 = x to the radius of an hypersphere and the three x m coordinates to distances 
measured on the hypersphere surface; time, x°, will still be measured along a 
direction normal to all others. If the hypersphere radius is very large we will 
not be able to notice the curvature on everyday phenomena, in the same way as 
everyday displacements on Earth don't seem curved to us. The Universe as a 
whole will manifest the consequences of its hyperspherical symmetry; using the 
Earth as a 3-dimensional analogue of an hyperspherical Universe, although our 
everyday life is greatly unaffected by Earth's curvature the atmosphere senses 
this curvature and shows manifestations of it in winds and climate. What we 
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propose here is an exercise consisting of an arbitrary assignment between coor- 
dinates and physical entities; the validity of such exercise can only be judged by 
the predictions it allows and how well they conform with observations. 

Hyperspherical coordinates are characterized by one distance coordinate, x 
and three angles p, 0, <p; following the usual procedure we will associate with 
these coordinates the frame vectors {Oi,G p ,O0,a(p}. The position vector for 
one point in 5D space is quite simply 

x = to t + xo % . (8.1) 

In order to write an elementary displacement dx we must consider the rotation 
of frame vectors, but we don't need to think hard about it because we can extend 
what is known from ordinary spherical coordinates. 

cbc = c7odf + a4dT + T<7pdp + Tsinpc7ed0 + Tsinpsin0a ( pd(p. (8.2) 

Just as before, we consider only null displacements to obtain time intervals; 

dt 2 = dx 2 + x 2 [dp 2 + sin 2 p (d0 2 + sin 2 0d<p 2 )] . (8.3) 

The velocity vector, v = x — Oq, can be immediately obtained from the displace- 
ment vector dividing by dt 

v = Ooi + xOpp + xsmpOo6 + xsmpsm60(p(p. (8.4) 

Geodesies of flat space are naturally straight lines, no matter which coordi- 
nate system we use, however it is useful to derive geodesic equations from a 
Lagrangian of the form (15.31) : in hyperspherical coordinates the Lagrangian be- 
comes 

2L = v 2 = x 2 + x 2 [p 2 + sin 2 p (e 2 + sin 2 6<p 2 )] . (8.5) 
Because de Lagrangian is independent of (p we can establish a conserved quantity 

J v = T 2 sin 2 p sin 2 0(p. (8.6) 

It may seem strange that any physically meaningful relation can be derived 
from the simple coordinate assignment that we have made, that is, proper time 
is associated with hypersphere radius and the three usual space coordinates are 
assigned to distances on the hypersphere radius. This unexpected fact results 
from the possibility offered by hyperspherical coordinates to explore a symmetry 
in the Universe that becomes hidden when we use Cartesian coordinates. In the 
real world we measure distances between objects, namely cosmological objects, 
rather than angles; we have therefore to define a distance coordinate, which is 
obviously r = xp. It does not matter where in the Universe we place the origin 
for r and we find it convenient to place ourselves on the origin. 
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Radial velocities r measure movement in a radial direction from our observa- 
tion point; we are particularly interested in this type of movement in order to find 
a link to the Hubble relation. Applying the chain rule and then replacing p 

t 

r = px + px = -r + px. (8.7) 

X 

We expect objects that have not suffered any interaction to move along o T ; from 
(18.41) we see that this implies p = 6 = <p = and then x becomes unity. Replacing 
in the equation above and rearranging 

t= 1 -. (8.8) 
r X 

What this equation tells us is exactly what is expressed by the Hubble relation. 
The value of x can be taken as constant for any given observation because the 
distance information is carried by photons and these preserve proper time, as 
we have seen in our discussion about electromagnetic waves. 1 The first member 
of the equation is the definition of the Hubble parameter and we can then write 
H = 1 jx. In this way we find the physical meaning of coordinate x as being the 
Universe's age. 

Underlying the present discussion there is an assumption a preferred frame 
where stillness means moving along a T ; there is no question of equivalent in- 
ertial frames here. This preferred frame is obviously attached to the observable 
still objects in the Universe which are galaxy clusters, as much as we can tell. 
This is far from the orthodox point of view, because galaxy clusters are seen as 
moving relative to each other and so cannot possible define a fixed frame. But 
in our formulation still objects move in straight lines along the proper time di- 
rection and keep their angular separations constant; this is naturally perceived 
as increasing mutual distances. If there is any relation between our formulation 
and an ether it must be found in the fact that movement has an absolute meaning, 
so it is defined relative to something that is fixed; we call the fixed reference a 
preferred frame while other authors call it ether. 

How does the use of hyperspherical coordinates affect dynamics in our labo- 
ratory experiments? We would like to know if these coordinates need only be 
considered in problems of cosmological scale or, on the contrary, there are im- 
plications for everyday experiments. The answer implies rewriting ( 18.21) with 
distance rather than angle coordinates; replacing p, 

dx = Oodt+ (o4--Op)dx + o p dr + r(o0d9 + sm6O(pd(p). (8.9) 



In order to preserve proper time photons must travel on the hypersphere surface and thus don't 
follow geodesies. 
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Evaluating time intervals from the null displacement condition, as before 



dt 1 



dT 2 - 2- drdr + dr 2 + r 2 (d0 2 + sin 2 0d<p 2 ) . (8.10) 

X 



This would be a version of ( 13.81) in spherical coordinates, were it not for the extra 
terms with powers of r/r in the second member. The coefficient rjx implies a 
comparison between the distance from the object to the observer and the size 
of the Universe; remember that x is both time and distance in non-dimensional 
units. We can say that ordinary special relativity will apply for objects which 
are near us, but distant objects will show in their movement an effect of the 
Universe's hyperspherical nature. 

With Eqs. (16.161) we have established the refractive indices n r and to account 
for the dynamics near a massive sphere using Cartesian coordinates; since this is 
frequently applied on a cosmological scale, we must find out how the dynamics 
is modified by the use of hyperspherical coordinates. Using the refractive indices 
and hyperspherical coordinates, noting that n r = n|, Eq. ( 14.71) becomes 



dt z = n\dx L + njx 2 dp 



(8.11) 



Dividing both members by dt 2 and reversing the equation 

njx 2 + n 4 4 x 2 p 2 =l; 



(8.12) 



and replacing xp by r — rx/x 



n\x 2j rn\ 



r 2 + ( - ) r L -2xf- 
XI X 



(8.13) 



Dividing both members by n\r 2 and rearranging results in the equation 



1 



i\ n 2 J 



1 



X \ xr 
- +2—. 

XI xr 



(8.14) 



As a further step we take the refractive index coefficients from Schwarzschild's 
metric (14.211) or those of from the exponential metric (16.161) and expand the sec- 
ond member in series of M/r taking only the two first terms. 



1-t 2 (2i 2 -4)M fx\ 



+ 



+ 2—. 
X J xr 



(8.15) 



The previous equation applies to bodies moving radially under the influence of 
mass M located at the origin which is, remember, the observer's position. For 
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comparison we derive the corresponding equation in Cartesian coordinates; start- 
ing with (18.121) it is now 

n 2 4 i 2 +njf 2 = l; (8.16) 
dividing by n\r 2 and rearranging 

1 x 2 \l ^1-t 2 (2i 2 -4)M 

-X 2 I ~2 2 1 3 ■ W.U) 

n 4 n %) r r r 

If we want to apply these equations to cosmology it is easiest to follow the 
approach of Newtonian cosmology, which produces basically the same results 
as the relativistic approach but presumes that the observer is at the centre of the 
Universe em. In order to adopt a relativistic approach we need equations 
that replace Einstein's in 4DO. A set of such was proposed above Eq. (16.41) but 
their application in cosmology has not yet been tested, so we will have to defer 
this more correct approach to future work. The strategy we will follow here is to 
consider a general object at distance r from the observer, moving away from the 
latter under the gravitational influence of the mass included in a sphere of radius 
r. If we designate by \x the average mass density in the Universe, then mass M 
in (18.151) is 47z:/ir 3 /3; this will have to be considered further down. 

Friedman equation governs standard cosmology and can be derived both from 
Newtonian and relativistic dynamics, with different consequences in terms of 
the overall size of the Universe and the observer's privileged position. From the 
cited references we write Friedman equation as 



rV 87T A 
r) =T" + 3 r 



— jU + T --; (8.18) 



with A a cosmological constant and k the curvature constant; the gravitational 
constant was not included because it is unity in non-dimensional units and the 
equation is written in real, not comoving, coordinates. In order to compare (18.151) 
with Friedman equation there is a problem with the last term because the Hubble 
parameter r/r does not appear isolated in the first member; we will find a way 
to circumvent the problem later on but first let us look at what (18.151) tells us 
when the mass density is zeroed. In this case = 1 and we find from (18.121) that 
t is unity, unless p is non-zero, for which we can find no reasonable explana- 
tion. Replacing and t with unity in (18.151) we find that r/r = 1/t, confirming 
what had already been found in (18.81) . Comparing with Friedman equation, this 
corresponds to a flat Universe with a critical mass density /x = fJL c ; it is imme- 
diately obvious that /i c = 3/(8ttt 2 ). Let us not overlook the importance of this 
conclusion because it completely removes the need for a critical density if the 
Universe is flat; remember this is one of the main reasons to invoke dark matter 
in standard cosmology. Notice also that this conclusion does not depend on a 
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privileged observer, because it is just a consequence of space symmetry and not 
of dynamics. 

Let us now see what happens when we consider a small mass density; here 
we are talking about matter that is observed or measured in some way but not 
postulated matter. The matter density that we will consider is of the order of 
1% of the presently accepted value. It is therefore just a perturbation of the flat 
solution that we described above and the fact that we are presuming a privileged 
observer has to be taken just for this perturbation. The first thing we note when 
we consider matter density is that t < 1, because there is now a component of the 
velocity vector along o p . Ideally we should solve the Euler- Lagrange equations 
resulting from (18.121) in order to find t and p but this is a difficult process and 
we shall carry on with just a qualitative discussion. Considering that we are dis- 
cussing a perturbation it is legitimate to make r/r « xjt and the two last terms 
in the second member of (18.151) can be combined into one single term (t/t) 2 , 
the same as we encountered for the flat solution, albeit with a numerator slightly 
smaller than unity. The first term has now become slightly positive and we can 
see from Friedman equation that this corresponds to a negative curvature con- 
stant, k, and to an open Universe. Lastly the second term includes the mass M 
of a sphere with radius r and can be simplified to 8^/i(t 2 — 2) /3; this has the 
effect of a negative cosmological constant; the combined effect of the two terms 
is expected to close the Universe W23L l24ll . The previous discussion was done 
in qualitative terms, making use of several approximations, for which reason we 
must question some of the findings and expect that after more detailed examina- 
tion they may not be quite as anticipated; in particular there is concern about the 
refractive indices used, which were derived in Cartesian coordinates both by the 
author and those that preceded him in using an exponential metric; it may happen 
that the transposition to hyperspherical coordinates has not been properly made, 
with consequences in the perturbative analysis that was superimposed on the flat 
solution. The latter, however, is totally independent of such concerns and al- 
lows us to state that the assumption of hyperspherical symmetry for the Universe 
dispenses with dark matter in accounting for the gross of observed expansion. 

Dark matter is also called in cosmolog y to account for the extremely high 
rotation velocities found in spiral galaxies 1E5I E^l and we will now take a brief 
look at how hyperspherical symmetry can help explain this phenomenon. Galaxy 
dynamics is an extremely complex subject, which we do not intend to explore 
here due to lack of space but most of all due to lack of author's competence to 
approach it with any rigour; we will just have a very brief outlook at the equation 
for flat orbits, to notice that an effect similar to the familiar Coriollis effect on 
Earth can arise in an expanding hyperspherical Universe and this could explain 
most of the observed velocities on the periphery of galaxies. Let us recall (18.91) . 
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divide by dt and invoke null displacement to obtain the velocity 

v= i + <Jpf+r(o e d + smO<J(p<p). (8.19) 

If orbits are flat we can make 6 = n/2 and the equation simplifies to 

v = ia 4 + (r- Op +r<pO(p. (8.20) 

Suppose now that something in the galaxy is pushing outwards slightly, so that 
the parenthesis is zero; this happens if r/r = i/x and can be caused by a pres- 
sure gradient, for instance. The result is that (18.201) now accepts solutions with 
constant r<p, which is exactly what is observed in many cases; swirls will be 
maintained by a radial expansion rate which exactly matches the quotient t/r. 
In any practical situation t will be very near unity and the quotient will be virtu- 
ally equal to the Hubble parameter; thus the expansion rate for sustained rotation 
is r/r « H. If applied to our neighbour galaxy Andromeda, with a radial extent 
of 30 kpc, using the Hubble parameter value of 81 km s _1 /Mpc, the expansion 
velocity is about 2.43 km s _1 ; this is to be compared with the orbital velocity 
of near 300 kms~ ! and probably within the error margins. An expansion of this 
sort could be present in many galaxies and go undetected because it needs only 
be of the order of 1 % the orbital velocity. 



9 Symmetries of G 4 1 algebra 

In this algebra it is possible to find a maximum of four mutually annihilating 
idempotents, which generate with an additive group of order 16; for a demon- 



stration see Lounesto [ 14, section 17.5]. Those idempotents can be generated by 
a choice of two commuting basis elements which square to unity; for the moment 
we will use 0023 an d O014. The set of 4 idempotents is then given by 



, (1+0-023) (l + o-ou) , O+O023XI-O014) 

h = ^ , h = 4 , 

, (1 -0023XI -0014) , (I-0023XI + 0014) 

h = 2 ' h = a • 



(9.1) 



Using the matrices of Sec.|3]to make matrix replacements of O023 and O014 one 
can find matrix equivalents to these idempotents; those are matrices which have 
only one non-zero element, located on the diagonal and with unit value. 

SU (3) symmetry can now be demonstrated by construction of the 8 generators 
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Ai = Ob2 (/l +/2) 



A 2 
A3 

A4 

A 5 



Ob3(/l+/ 2 ) 



(73 + O02 

2 ' 

-a 2 + oq3 



/1-/2 



C7Q14 — <7i234 



-CTl(/2 + /3) 
-04(/ 2 + /3) 
O(012(/l 4-/3) 



A7 = - O024 (/l 4-/3) 



-<7i -C7 4 



-C74 + C701 



O012 4- O034 



Ob 13 — (7024 



A 8 



f\+fl — 2/3 2ao23 + C7Q14 4" C71234 



x/3 



2^ 



These have the following matrix equivalents 



Xa = 



/o 1 0\ 

10 



\0 0/ 

/o 0\ 
10 
10 

\o oy 

/0 -j 


j 

\0 



/o 
j 





A 7 = 



A 5 

o\ 






/o 



j 
\0 



-j 0\ 





0/ 

0\ 

-j 






A,= 





0/ 



A 8 =(l/v / 3' 



A 6 = 

/ 1 

1 

\0 



(\ 0\ 
0-100 



\0 0/ 

/o 1 0\ 



10 

\0 0/ 

0\ 










0/ 



(9.2) 



(9.3) 



which reproduce Gell-Mann matrices in the upper-left 3*3 corner H15LI27L 12811 . 
Since the algebra is isomorphic to complex 4*4 matrix algebra, one expects to 
find higher order symmetries; Greiner and Miiller [27] show how one can add 7 
additional generators to those of SU (3) in order to obtain SU (4) and the same 
procedure can be adopted in geometric algebra. We then define the following 
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additional SU (4) generators 



A9 

A10 
An 
A12 

A13 
A14 

A15 



^l(/l+/4) = 
04(/i+/ 4 ) = 
— C012(/2 + /4 
^024(/2 + /4) : 

o- 3 (/ 3 +/ 4 ) = 

<72(/3+/4) = 



^1 ~ <7Q4 
2 

04 + obi 



-0012 — 0034 



_ 0013 + <7Q24 

~ 2 : 

03 ~ <7Q2 

2 ' 

02 + CT03 



(9.4) 



/l + H + h - 3/4 O023 - C014 - CTi234 



V6 



V6 



Once again, making the replacements with Eq. (13.91 ) produces the matrix equiv- 
alent generators 



Aq = 



A 



12 



/0 1\ 





\1 0/ 

/o 0\ 
-j 


V° j °/ 



A10 = 



An 



A 



13 



45 



/0 -j\ 



\j 00 0/ 

/o 0\ 



1 

\0 1 0/ 

/l \ 

10 

10 

\0 -3/ 



A 



14 



/o 0\ 

1 



\0 1 0/ 

/o 0\ 

-j 

Vo j 0/ 



1/V6 



(9.5) 



The standard model involves the consideration of two independent SU (3) 
groups, one for colour and the other one for isospin and strangeness; if gen- 
erators Ai to As apply to one of the SU (3) groups we can produce the generators 
of the second group by resorting to the basis elements 03 and Oo4- The new set 
of 4 idempotents is then given by 



_ (l+<T 3 )(l+go 4 ) _ (l+g 3 )(l-Ob4) 

Jl 4 , J2 4 ' 

(I-CT3XI-CT04) r (l-C7 3 )(l+O04) 

h = -. j J4 = • 



(9.6) 
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Again a set of SU (3) generators can be constructed following a procedure similar 
to the previous one 



™ ^ /r , m °02 + O023 

C01 + C013 



a2 = o- i(/i+/ 2 

«3=/l-/2 



2 

004 — O034 



a 4 = o- 2 (/2 + /3) = 

«5 = -Oi(/2+/3) 
06 = 04(/l+/3) = 

a7 = o- i 2 (/i+/ 3 ) 
a 8 



2 

(72 + <7Q24 
2 

-CTl ~ gQ14 (9 - 7) 

2 

O4 ~ <7Q3 

2 ' 

0012 + 01234 



2 

fl+fl- 2/3 2(73 + O04 + 0034 



V3 2V3 

This new 5t/ (3) group is necessarily independent from the first one because its 
matrix representation involves matrices with all non-zero rows/columns, while 
the group generated by X\ to X% uses matrices with zero fourth row/column. In 
the following section we will discuss which of the two groups should be associ- 
ated with colour. 

At the end of Sec.|3]we used one particular idempotent to split the wavefunc- 
tion into left and right spinors and here we discuss how the different idempotents 
are related to the symmetries discussed above, suggesting a relation between 
idempotents and the different elementary particles. We have already established 
that each set of 4 idempotents is generated by a pair of commuting unitary basis 
elements. Let any two such basis elements be denoted as h\ and h,2\ then the 
product /?3 = h\h,2 is itself a third commuting basis element. For consistence we 
choose, as before, 

h\ = O023, h 2 = O014; (9.8) 

to get 

h = <?1234, (9.9) 

which commutes with the other two as can be easily verified. The result of this 
exercise is the existence of triads of commuting unitary basis elements but no 
tetrads of such elements. We are led to state that a general unitary element is a 
linear combination of unity and the three elements of one triad 

/? = ao + ai/?i +a2^2 + ^3^3- (9.10) 
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Since h is unitary and the three h m commute we can write 

h 2 = [(ao) 2 + (ai) 2 + (a2) 2 + (a3) 2 ]+2(aoai-a 2 a3)hi + 
+ 2(ao«2 — a\a?,)h2 + 2(ao«3 — a\a2)h^ = 1 

The only form this equation can be verified is if the term in square brackets is 
unity while all the others are zero. We then get a set of four simultaneous equa- 
tions with a total of sixteen solutions, as follows: 8 solutions with one of the 
equal to ±1 and all the others zero, 6 solutions with two of the equal to —1/2 
and the other two equal to 1/2 and 2 solutions with all the simultaneously 
±1/2. The coefficients play the role of quantum numbers which determine 
the particular idempotent that goes into Eq. (13.171) : these unusual quantum num- 
bers are expressed in terms of the SU (4) generators A3, As and A15 in Table[T]in 
order to highlight the symmetries. We don't propose here any direct relationship 



Table 1: Coefficients for the various unitary elements. 
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O1234 


a 3 


A 8 


A15 


(«o) 


(fli) 
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1/V3 
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-1/2 
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-1/2 


-1/2 


1/2 


-1 


-1/V3 
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-1/V3 


-1/V6 
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1/2 


-1/2 
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0/2 


1/2 


1/2 


1/2 


1/2 





2/V3 


-1/V6 


-1/2 


-1/2 


-1/2 


-1/2 





-2/V3 


1/V6 



between the various idempotents and the known elementary particles, although 
the fact that the standard model gauge symmetry group is found as direct conse- 
quence of the monogenic condition which itself generates the Dirac equation is 
rather intriguing. 
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10 Conclusion and future work 



Monogenic functions applied in the algebra of 5-dimensional spacetime have 
been shown to originate laws of fundamental physics in such diverse areas as 
relativistic dynamics, quantum mechanics and electromagnetism, with possible, 
still unclear, consequences for cosmology and particle physics. To say that those 
functions provide us with a theory of everything is certainly unwarranted at this 
stage but it is clear that there is a case for much greater effort being invested in 
their study. 

There are unanswered questions in the present work. For instance, how can 
we avoid an ad hoc definition of inertial mass or what is the true relation between 
the symmetries generated by monogenic functions and elementary particles? In 
spite of its various loose ends, the formalism is perfectly capable of unifying 
relativistic dynamics, quantum mechanics and electromagnetism, which in itself 
is no small achievement. Certain developments seem relatively straightforward 
but they must be made, even if no knew predictions are expected. Applying 
monogenic functions to the Hydrogen atom should not be difficult because the 
form of the Dirac equation we arrived at is perfectly equivalent to the standard 
one; one should then find the same solutions but in a GA formalism. In the 
same line one could try to solve the equation for a central gravitational potential, 
being certain to find quantum states. It is not clear how important these could be 
in planetary mechanics or galaxy dynamics. 

Gravitational waves are predicted by the monogenic function formalism as we 
pointed out but did not investigate. How important are they and what chance is 
there of them being detected by experiment? We don't know the answer and we 
don't know what difficulties lie on the path of those who try to solve the equa- 
tions; this is an open area. The sources tensor must be clearly understood and 
directly related to geometry; at the moment all densities, mass, electromagnetic 
energy, etc. must be inserted in the equations but one would expect that a per- 
fect theory would produce such densities out of nothing. In previous papers we 
suggested that a recursive, non-linear, equation could be the answer to the prob- 
lem but the concept has not yet been formalized and there are no clear ideas for 
achieving such goal. 

In conclusion, the present work opens the gate of a path that will possibly lead 
to an entirely new formulation and understanding of physics but this path is very 
likely to have many hurdles to jump and several dead ends to avoid. 

A Indexing conventions 

In this section we establish the indexing conventions used in the paper. We deal 
with 5-dimensional space but we are also interested in two of its 4-dimensional 
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subspaces and one 3-dimensional subspace; ideally our choice of indices should 
clearly identify their ranges in order to avoid the need to specify the latter in 
every equation. The diagram in Fig. [T] shows the index naming convention used 
in this paper; Einstein's summation convention will be adopted as well as the 




Figure 1: Indices in the range {0,4} will be denoted with Greek letters a,/3, y. 

Indices in the range {0,3} will also receive Greek letters but chosen 
from /!, v,£. For indices in the range {1,4} we will use Latin letters 
i, j, k and finally for indices in the range {1,3} we will use also Latin 
letters chosen from m, n, o. 

compact notation for partial derivatives d a = d / dx a . 

B Time derivative of a 4-dimensional vector 

If there is a refractive index the wave displacement vector can be written as 

dx = g a dx a = nP a o p dx a . (B.l) 

Because this vector is nilpotent, by virtue of Eq. (14.61) . the five coordinates are not 
independent and we can divide both members by dx° = dt defining the nilpotent 
vector 

x = go + gix 1 = n a Q o a +nP jOpx' . (B.2) 

Suppose we have a 5D vector a = o a a a and we want to find its time derivative 
along a path parameterized by t, that is all the x' are functions of t. We can write 

a = dpa a xPo a \ (B.3) 

where naturally x = 1. Remembering the definition of covariant derivative 
(12.251) and Eq. (IB .21) we can modify this equation to 

ci=xPg j3 -g l5 dpa (X Oa=x-(Da). (B.4) 

We have expressed vector a in terms of the orthonormed frame in order to avoid 
vector derivatives but the result must be independent of the chosen frame. 
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This procedure has an obvious dual, which we arrive at by defining 



x = g^+g 4 . 



(B.5) 



The proper time derivative of vector a is then 



a = x- (Da). 



(B.6) 
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